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We introduce a linear, canonical transformation of the fundamental single-mode field opera- 
tors a and a* that generalizes the linear Bogoliubov transformation familiar in the construction of 
the harmonic oscillator squeezed states. This generalization is obtained by adding to the linear 
transformation a nonlinear function of any of the fundamental quadrature operators X\ and X2, 
{— ^ making the original Bogoliubov transformation quadrature-dependent. Remarkably, the conditions 

■ of canonicity do not impose any constraint on the form of the nonlinear function, and lead to a 

' set of nontrivial algebraic relations between the c-number coefficients of the transformation. We 

examine in detail the structure and the properties of the new quantum states defined as eigenvec- 
tors of the transformed annihilation operator b. These eigenvectors define a class of multiphoton 
squeezed states. The structure of the uncertainty products and of the quasiprobability distributions 
in phase space shows that besides coherence properties, these states exhibit a squeezing and a de- 
formation (cooling) of the phase-space trajectories, both of which strongly depend on the form of 
the nonlinear function. The presence of the extra nonlinear term in the phase of the wave functions 
has also relevant consequences on photon statistics and correlation properties. The non quadratic 
structure of the associated Hamiltonians suggests that these states be generated in connection with 
^ , multiphoton processes in media with higher nonlinearities. 

^ ! PACS. 03.65.Sq, 42.50.Dv. 
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I. INTRODUCTION 



The generalization of the squeezed states of the harmonic oscillator to different systems and the definition of new 
pi ' nonclassical states of light have been actively pursued in recent years | pj, , Q| , B, f| , Q] . The possibility of 
experimental realization has also been discussed in some instances ||] , Q , |1C| ] . However, these proposed generalizations 
cannot be easily and straightforwardly connected to canonical and Hamiltonian structures. Generalizations providing 
& \ a Hamiltonian structure and at the same time preserving the canonical commutation relations would thus be extremely 
desirable not only from a mathematical and conceptual point of view, but expecially for possible connections to new 
physical systems and experimental realizations. 

In this paper we present a natural generalization of the Bogoliubov transformation originally introduced by Yuen to 
define the two-photon coherent states of the harmonic oscillator [|ll | . This generalization is achieved by simply adding 
to the original transformation, which is linear in the fundamental mode variables a and , a nonlinear operator- 
valued function F(Xi) of any of the fundamental quadratures X\ and X2. The nonlinear function F is by definition 
hermitian, and thus defines a quadrature-dependent Bogoliubov transformation. The transformed modes b and W obey 
the canonical commutation relations [b, b^] under suitable constraints on the coefficients of the transformation. Such 
constraints are in fact broad enough to allow for enough freedom on the variation of the coefficients. The eigenvectors 
of the transformed annihilation operator b define the quadrature-dependent generalization of the squeezed states of 
the harmonic oscillator. Such states are thus multiphoton squeezed states associated to non quadratic Hamiltonian 
systems. They preserve some coherence properties, such as the classical motion of the wave-packet center, but in 
addition they exhibit nonclassical features such as amplitude squeezing, interference and deformation in phase space, 
super-Poissonian modifications of the field statistics and of the second-order correlation functions. To avoid possible 
sources of confusion, we emphasize that the multiphoton squeezed states defined in the present work should not be 
confused with the nonlinear coherent states defined in the literature [jl0| as the eigenstates of nonlinear operator 
valued functions of the number operator. In fact, both the physical origin and the mathematical properties of these 
two classes of states arc rather different; they are both associated to multiphoton down conversion processes, but the 
multiphoton squeezed states are generated in multiphoton down conversion processes in optical parametric amplifiers 
of higher x nonlincarity, while the nonlinear coherent states are produced in multiphoton down conversion processes 
in the quantized motion of a trapped atom (the counterpart of this process has also been recently studied [ p"2] ) . 
The plan of the paper is the following: the quadrature-dependent Bogoliubov transformation is introduced in Sec. 
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II, and the conditions of canonicity are discussed. In Sec. Ill we introduce and discuss in some detail the Hamiltonian 
structure induced by the transformation; in particular, we provide the explicit form of the Hamiltonian expressed in 
terms of the original mode variables for the lowest nonlinear function of the amplitude F(Xi) = X\. The properties of 
the uncertainty products are studied in Sec. IV, where it is shown that the multiphoton squeezed states associated to 
the quadrature-dependent Bogoliubov transformation exhibit amplitude squeezing and are moreover asymptotically 
of minimum uncertainty. The eigenvalue equation defining the multiphoton squeezed states are explicitly solved in 
the quadrature representation in Sec. V, and the coherent motion of the wave-packet center is obtained. Results 
for the photon number distributions and second-order correlation functions are presented in Sec. VI, while in Sec. 
VII we perform the phase-space analysis of the multiphoton squeezed states in terms of the Wigner quasiprobability 
distributions and of the classical parametric trajectories. Comments and conclusions follow in Sec. VIII. 



II. THE QUADRATURE-DEPENDENT BOGOLIUBOV TRANSFORMATION 

Let us consider a single-mode bosonic quantum system. The fundamental adimensional variables are the bosonic 
annihilation and creation operators a and a\ with canonical commutation relation [a, a^] = 1. 
We introduce a new set of fundamental variables b and b^ by the linear transformation 

b = M a + i/at + iF{X 1 )+ rjG(X 2 ) , (1) 

and its adjoint 

6 1 = n*cJ + v*a + 1 *F{X 1 ) + n*G{X 2 ) . (2) 

In the above, [i, v, and 7 are complex numbers, while F — and G — are operator valued nonlinear functions 
of the fundamental quadrature operators 

q + qt a - a t . . 

*i = -2-, ■ (3) 

The transformation is obviously a linear mapping in the space of operators. On the other hand, it is in general a 
nonlinear, quadrature-dependent extension of the Bogoliubov transformation in the sense that the latter is recovered 
either if 7 = 77 = 0, or if F and G contain only linear powers of the quadratures. 

To begin with, we will study, for simplicity, only the quadrature-dependent Bogoliubov transformation generated 
by the first quadrature X\ alone, i.e. we put 77 = in Eqs. (|l|) and (||) and thus consider the transformation 

b = /ia + z/a f +7F(Xi) , (4) 

together with its adjoint 

=p*a t + v*a + 1 *F {X x ) . (5) 

In fact, all the main features of the quantum states generated by the quadrature-dependent Bogoliubov transfor- 
mations are already contained in this instance. The extension to transformations generated by G{X 2 ) alone, and by 
both F(X\) and G(X 2 ) will be briefly discussed in the conclusions. 

i) Canonical structure. We first analyze under what conditions the quadrature-dependent Bogoliubov transforma- 
tion is canonical. With a little algebra, we have that [b, tf\ = 1 if 

M 2 -M a = i. 



Re(/J7* - = , (6) 

where Re(z) denotes the real part of the complex number z. It is a remarkable feature that the conditions of canonicity 
for the quadrature-dependent Bogoliubov transformation are relations for the c-number coefficients of the transfor- 
mation only, and do not involve the operator- valued function F{X{). The same is true also when considering the 
quadrature-dependent Bogoliubov transformations generated by G(X 2 ) alone, or by both F{X{) and G(X 2 ). There- 
fore, the form of the operator-valued functions F(X{) and G(X 2 ) is not constrained by requiring the transformation 
to be canonical. 

The relations are such to give enough freedom in the choice of the coefficients of the transformation. This 
would not be true if we had attempted to generalize the linear Bogoliubov transformation introducing nonlinear 
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operator-valued functions A(a) and/or -B(a^) of the fundamental mode variables rather than of the quadratures. It 
is easy to verify that in this case the conditions of canonicity lead to such strict constraints on the coefficients of the 
transformation that, in practice, only the linear case 7 = is attainable. Thus it is the requirement of canonicity that 
selects the quadratures as the arguments of the nonlinear operator-valued functions in the canonical transformations. 

We see that the first of relations (Q) is exactly the same as in the linear case (7 = 0), but it is now constrained to 
be compatible with an additional relation, the second of Eqs. ([|), which connects 7, the cofficient of the nonlinear 
part of the transformation with the coefficients \i and v of the linear part. 

ii) Eigenstate of b: the multiphoton squeezed state. We now introduce, in strict analogy with the linear case, the 
cigenstatc | "J) p of the transformed annihilation operator b: 

b\^) = m)0, (7) 

where (3 is a c-number. The eigenvalue equation (]?]) defines a new class of quantum states, the multphoton squeezed 
states, which are a direct generalization of the two-photon coherent states originally introduced by Yuen via the 
linear Bogoliubov transformation. In principle, depending 011 the choice of F(Xi), the Bogoliubov nonlinear squeezed 
states include all different generalized many-photon coherent states, from four-photon coherent states on, as shall be 
illustrated below. 



III. HAMILTONIAN STRUCTURE 



As already mentioned, an important feature of the new quantum states defined by Eq. ^ is that, thanks to 
the canonical quadrature-dependent Bogoliubov transformation (U)-(|5|), they can be associated to a non quadratic 
Hamiltonian structure in terms of the original mode variables a and . In fact, the quadratic Hamiltonian in the 
transformed variables 



H = tfb (8) 



reads, in terms of the original mode variables, 



H = (|/i| 2 + M Va + M*^ f2 + ^*a 2 + M 2 + \ 

+ /i* 7 a t J F 1 (X 1 ) + ^*F{X 1 )a + v<y*F(Xx)a,i + v*jaF(Xx) + |7| 2 F 2 (Xi) . (9) 

The first line of the Hamiltonian (^|) contains the original Yuen squeezing Hamiltonian, including the two-photon 
contributions a' 2 and a 2 . The remaining terms appearing in the second line introduce anharmonicities, as they involve 
higher powers of a and , whose degree depends on the form of F{X\). For instance, if we choose the lowest nonlinear 
power, i.e. F(Xi) = X 2 , additional linear, quadratic, cubic, and quartic terms in a and are generated. We are 
thus dealing with the effective description of a single-mode quantum system undergoing multiphoton down conversion 
processes (parametric down conversion processes of higher order). The full Hamiltonian in the case F(X\) 

H=^- (at 4 + a 4 ) + | 7 | 2 (a^a + a V) + i| 7 | VV 

+ \ (M*7 + ^7*) a t3 + i (/i7* + a 3 + Re [(fi + ^7*] (a t2 a + a) a 2 ) 
bd! . „*„\ „ta . ftL _u ,,,A „2 -u fl„l 2 -l. M 2 . OUR „t, 



U±- + M Vj a t2 + IU±- + a 2 + (| M | 2 + M 2 + 2| 7 |-j a<a 

+ l - (m* 7 + 3^7* + 2^*7) a) + l - ( M7 * + 2^7* + 3^* 7 ) « + ^ + W\ 2 + \ ■ (10) 

In general, choosing F(X\) — X\ L with n > 2 will generate a 2n-photon Hamiltonian. However, if we think of 7 as 
the strength of the effective coupling to the optical medium with higher nonlinearity, we sec from Eqs. (§)-© that 
in the typical weak-coupling situations that are realistically foreseeable the dominant terms will be those associated 
to the lowest powers of 7. One can of course consider trascendental forms for F(Xi), such as, e.g., F = sinXi; also 
in these cases we expect that in the weak-coupling regimes the lowest terms in a power series expansion of F will 
contain the dominant contribution from the nonlinear part of the Hamiltonian. 
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It is instructive to write the Hamiltonian (||) in terms of the quadratures X\ and X2, as this clarifies the physical 
effects associated to the different terms. We adopt, due to canonicity, the standard parametrization /1 = coshre 1 *' 1 , 
v = sinhre 1 ^ 2 . We have (letting 7 = |7|e iA ): 



H = [e 2r (l + cos (fa - fa)) + e- 2r (l 



fa))} ^ + [e 2r (l - cos (fa - fa)) + e- 2r (l + cos (fa - fa))] ^ 



sinh 2r sin I 



i){X 1 ,X 2 } + 2V2\-/\coshrcos(fa-8)X 1 -F(X 1 ) 



hi 

V2 



[coshrsin(0i - 5) - sinh r sin (fa - 5)]{F(X 1 ),X 2 } + |7| 2 F 2 (Xi) 



(11) 



where { , } denotes the anticommutator. We see from Eq. (|flj) that the Hamiltonian contains the standard terms 
associated to squeezing; however the additional terms generated by -fF(Xi) clearly modify the distribution in phase 
space with respect to the standard two-photon squeezing, and this in turn affects both the structure of the states 
and the observable properties, as will be seen below. In fact, the quadrature-dependent Bogoliubov transformation 
really amounts to a drastic deformation of the second quadrature X 2l i.e. the momentum in particle language. This 
can be most easily seen by considering the particular case fa- — fa — 4>- I n such a situation, due to the canonical 
constraints (||) the phase 4> of the coefficients (i and v and the phase 5 of the coefficient 7 must obey the relation 
<f) = S + ((2k + l)-7r)/2, and the transformed Hamiltonian reduces to 



(X 2 + 2\ 1 \e r F(X 1 )) i 



(12) 



The above expression shows immediately that in the limit 7-^0 the Hamiltonian (y_2|) reduces exactly to the 
standard squeezing Hamiltonian. We thus see that the canonical quadrature-dependent Bogoliubov transformation 
is an effective description of the coupling of the system to an external medium with the following overall effects: 
a squeezing transformation in both quadratures and a shift in the second quadrature (the canonical momentum in 
particle picture, the field phase in optical picture). The latter is enforced by the nonlinear function F(Xi), and for 
F = X 1 it reduces to the well-known Yuen shift in the momentum. Of course, the shift F(Xi) is not and must not 
be interpreted as an external electromagnetic vector potential. Instead, it summarizes the effect of the interaction 
with the external medium. In turn, the shift affects the squeezing properties according to the specific form of F(Xi) 
and modifies the structure of the stable, closed trajectories in phase space. We thus speak of quadrature-dependent 
squeezing and of quadrature-dependent shift and distorsion (deformation) of phase-space trajectories and phase-space 
distributions. However, before turning to a detailed analysis of the statistical properties (photon number distribution 
and correlation functions) and of the phase-space distributions we first discuss the structure of the uncertainty 
products and of the eigenstates of the transformation. 

IV. UNCERTAINTY PRODUCTS AND ASYMPTOTIC MINIMUM UNCERTAINTY 



Let us consider the quantities A 2 Xi = (Xf) — (Xi) 2 (i = 1,2). We want to evaluate the uncertainty product 
A 2 Xi A 2 X 2 in the multiphoton squeezed state \^)/3 solution of the eigenvalue equation (^). In terms of the transformed 
modes b and the fundamental quadratures read 



where 



X, 



(!) 



X 1 = \[( l i-v)% + ( l i-u)b^] , 



Xo=X, 



CO 



X 



(nl) 



[{fi + up-ift + uyb] , 



(13) 



Im(^7* — v*"f)F 



Gu-i/)*& + 0u-i/)&t 



(14) 



where Im(z) stands for the imaginary part of the complex number z. The expression for the first quadrature X\ 
is identical to the one obtained via the linear Bogoliubov transformation: this is obviously due to the fact that 
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the nonlinear part of the canonical quadrature dependent Bogoliubov transformation depends only on X\ . What is 
changed with respect to the linear case is evidently the expression for the second quadrature X- We can see that it 
is made of two terms: the first term X% is identical to that obtained through the linear Bogoliubov transformation, 

while the second term X^™^ is the contribution due to the nonlinear part of the transformation, where the argument 
X\ of the nonlinear operator-valued function F has been re-expressed using the the first of Eqs. (|l3|). 
Direct evaluation of the uncertainties in the state yields for the first quadrature 

A 2 X 1 = i| M -H 2 , (15) 

which is obviously still the same expression originally obtained by Yuen in the linear case. For the second quadrature 
one obtains 

A 2 X 2 = A 2 xf + A 2 X ( 2 nl) + ({x { 2 l \x { 2 rd) }) - 2{X { 2 l) ){X ( 2 nl) ) . (16) 

We see that the uncertainty in the second quadrature is substantially modified with respect to the linear case: the 
canonical quadrature-dependent Bogoliubov transformation introduces additional terms; in particular, the presence of 
the anticommutator expresses the existence of correlations between the linear and the nonlinear part of X . Obviously, 
the term A 2 X ' acquires the same value as in the linear case: 

A 2 X« = i| M + H 2 , (17) 

while, introducing (jx^X^W = ( jx^°, X^ n0 j) - 2(x! 2 J) )(X i 2 nl) ), the remaining terms in Eq. (|l|) read 



<{xf ,xf >}) T =Im( M 7* -^7)({xiV(*i)}>T , (18) 

where ({xj°, F(Xi)|) T = ({xf,F(X)}) - 2(Xj/ ) )(F(X)). At variance with the uncertainty @ in X$\ the 

remaining terms (|l^) that contribute to the overall uncertainty in X depend explicitely on the eigenstate \^)/3 through 
the eigenvalue (3. Collecting terms together, the uncertainty product reads 

A 2 XA 2 X = l|/i 2 - ^ 2 | 2 + ilm 2 ( M7 * - i/* 7 )| M - ^| 2 A 2 F(X) 
lb 4 

+ hm(^* */* 7 )|/x v\ 2 {{x$\ F(X 1 )}) T . (19) 

Setting 7 = in Eq. ( ]l9| ) and choosing (j>\ = (f>2 = 4> i n the expressions for the coefficients /i and v of the linear 
part of the transformation, the uncertainty product reduces to the standard Hciscnberg minimum with equal and 
opposite squeezing in the quadratures X and X- In the general case 7^0 the extra terms are in general non-zero. 
In the particular case [i and v both real or with equal phases, then A 2 X = (1/4) exp (— 2r), A 2 X 2 ^ = (1/4) exp (2r). 
Moreover, in this case, letting 7 = 71+172, the conditions of canonicity (^) imply Im(/i7* — v*~{) = —(72 exp (r))/ cos cj>. 

If we choose F(Xi) of the form F = Xf with the integer n > 2, then A 2 F ~ exp (-2nr) and (jx^, F(Xi)}) T ~ 
exp (— 2(n — l)r). Collecting terms together we finally have 

A 2 X A 2 X * 4 + 7l^ 2nr + l2e~ 2nr . (20) 
16 

We thus see that for large values of r and/or for small values of 72 the uncertainty product is close to the Heisenberg 
minimum, and we can speak of states of quasi minimum uncertainty. 
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V. SOLUTION OF THE EIGENVALUE EQUATION IN THE QUADRATURE REPRESENTATION 



i) Multiphoton squeezed vacuum. We first consider the equation for the vacuum state \^)o of the transformed mode 
b, i.e. for the multiphoton squeezed vacuum: 



6|*)o = 



(21) 



In the quadrature (Xi) representation the operator X\ acts as a multiplication by the c-number x, and the equation 
for the vacuum state reads 



(coshre^i) X + * X2 + (sinhre^ 2 



x — iXn . . 



*o(a0 = 



(22) 



where the wave function ^q(x) = (Xi\^)o is the multiphoton squeezed vacuum in the quadrature representation. We 
need to recall that the conditions of canonicity (||) impose 

cosh r cos (0i — S) — sinh r cos (02 — S) = . (23) 

Multiplying both sides of Eq. ( |22| ) by exp (— i<5), exploiting relation ( p^ ) and solving for ^o(x) we finally obtain 

^o(x) =^ 1/2 cxpf-^)exp 



( coshrcos(0i -S)x 2 + 2| 7 | / dx'F{x') 



(24) 



where Ao is the normalization factor, a = C_/C+ is the variance of the Gaussian density profile of the wave function 
and 



C+ = cosh r sin (0i — 5) + sinh r sin (02 — , 



cosh r sin I 



<5) — sinh r sin (02 — 5) 



(25) 



The structure of the squeezed vacuum is thus the following: the probability density is of Gaussian shape centered 
in (X\) = 0, while the function F(x) enters only in the phase of the wave function. The phase is made up by two 
contributions: a typical squeezing correlation term proportional to x 2 plus a term proportional to the integral of F(x). 
As the parameters of the transformation, and consequently the wave function, may in principle be time-dependent, one 
can easily derive the equations of motion for the center of the wave packet (Xi) (t) and for the squeezing parameter r(t). 
Considering the case 0i = 02 = one has = 5 ± ((2k + 1)tt)/2, i.e. cos (0 — 5) = 0, C+ = exp (r), C_ = exp (— r), 
a = exp (— 2r), and the wave function reduces to a Gaussian modulated by a phase factor which includes only the 
integral of the function F(x). 

ii) Multiphoton squeezed states. We turn now to the solution of the eigenvalue equation (J?]) that defines the 
multiphoton squeezed states for arbitrary complex values of the eigenvalue j3 = |/3|exp(i£). Introducing the wave 
function ^/3(x) = (Xl|\E') / 3 the eigenvalue equation (^) is easily solved in the general case along the same lines exposed 
in the case of the multiphoton squeezed vacuum and the solution reads 



*p(x)=Af- 1 ' 2 



exp 



(a: - x ) 2 






2(7 


exp 





l r cos ( 



5)x 2 + \/3\ cos (£ - S)x - 2| 7 | / dx'F(x') 



(26) 



where the normalization N^ 1 — l/y/ira, and Xo = (\(3\/C+) sin (£ — 5) is the center of the Gaussian probability 
profile. The phase of the wave function contains both the x-dependent and the x 2 -dependent terms typical of the 
squeezed coherent states of the harmonic oscillator, plus an additional anharmonic contribution coming from the 
F(x)-dependent term. In the particular case 0i — 02 = (i.e. 5 = ±((2fc + 1)tt)/2), introducing the complex 
number a = ct\ + icti we can choose /3 = fia + va* , which, for the harmonic oscillator, is the condition of equivalence 
between the Yuen two-photon coherent states and the standard squeezed coherent states obtained first by squeezing 
the vacuum and then by displacing it. We thus have, writing (3 = /?i + i(3-2, Pi = ct\ exp (r) for the real part, and 
fh = oil exp (— r) for the imaginary part. The wave function then reduces to 



^p(x) =A/' ( r 1 / 2 exp 



e 2r 

-— (x - \f2~a x 



exp 



\/2i 



a 2 x - e |7| 



dx'F(x') 



(27) 



G 



The state (27) describes a coherent squeezed dynamics: (Xi) = \pla\ is still equal to the mean in a coherent state 
^ a {x) of the harmonic oscillator, while the variance a = cxp (— 2r)/2 coincides with the dispersion A 2 Xi of the Yuen 
two-photon coherent state. The equation for the wave-packet center is easily derived, and one obtains 

d 2 (X 1 ) 



dt 2 



+ e"' (^(F(X 1 ))^ + \ 7 \(d t F(X 1 )) S j 



(28) 



which reduces to the equation of motion for the classical harmonic oscillator when neither 7 nor F(Xi) are explicitly 
time-dependent. Therefore the multiphoton squeezed states preserve some of the basic properties of the squeezed 
coherent states of the harmonic oscillator. However, the presence of the F-dependent term in the phase of the wave 
function has relevant effects on other physical quantities, in particular, as shall be seen in the next sections, on the 
photon statistics, the correlation properties and the structure of the quasiprobability distributions in phase space. We 
also note from the above discussion that in all cases F(x) enters only in the phase of the wave function: it is therefore 
always possible to cast the Hamiltonian (O) in the general form 

H = aXf + {bX l + cX 2 + 7^(^i)) 2 , (29) 

with coefficients a, b, c that depend on the parameters //, v, 7 of the canonical quadrature-dependent Bogoliubov 
transformation. 



VI. FIELD STATISTICS 



Let us consider the photon number distribution (PND) in a multiphoton squeezed state: P(n) — | (nl^)^ 2 . Due to 
the nonlinear nature of the function F(Xi) it is in general impossible to write a closed analytic expression for P(n), 
which can however be easily plotted numerically. To gain immediate insight on the effect of the F-dependent term 
that enters in the phase of the state $ ' p(x) we have confronted the PND of the standard Yuen two-photon coherent 
state with the PND of the multiphoton squeezed state with the lowest nonlinear behavior F(x) = x 2 . We compare at 
<pi = <p2 = 0, j3 real, and set the coupling 7 at intermediate small values. In Fig. 1 we have plotted the PND of the 
multiphoton squeezed state at 7 = 0.1, f3\ = 3 and squeezing parameter r = 0.8 versus the PND of the corresponding 
Yuen two-photon coherent state (7 = 0). 
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FIG. 1. The PND P(n) of the multiphoton squeezed state (^) (full line) plotted versus the PND of the corresponding 
two-photon coherent state (dotted line). For both PND's the eigenvalue /3i = 3 and the squeezing parameter r = 0.8. The 
PND of the multiphoton squeezed state has been evaluated for the lowest nonlinearity F(x) — x 2 and a coupling 7 = 0.1. 

We notice from Fig. 1 that the PND of the multiphoton squeezed state lowers the first maximum that appears 
in the PND of the corresponding two-photon coherent state at a very low number of photons. Correspondingly, 
at a higher number of photons the PND of the multiphoton squeezed state stays well above zero approximately up 
to n = 20 and well above the PND of the two-photon coherent state, which is practically zero already at n = 10. 
It is also to be noted that the presence of the nonlinearity tends to attenuate the oscillations in the PND, and we 
have tested that this effect becomes more pronounced for increasing values of the coupling 7. In fact, Fig. 1 shows 
that the minima of the oscillations in the PND of a multiphoton squeezed state stay always well above zero. The 
explanation of this phenomenon is that the additional cubic nonlinear term in the phase of the wave function yields 
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faster oscillations that enter in competition with the slower oscillations due to the remaining linear and quadratic 
squeezing terms entering in the phase. The net result is a decrease in the maxima and a corresponding increase in 
the minima of the oscillations. The effect is obviously more pronounced if we consider the PND of a multiphoton 
squeezed state with higher nonlinearity, as shown in Fig. 2 for the case F(X\) = Xf. 
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. 01 
.0075 

.005 
.0025 

20 40 60 80 100 

FIG. 2. The PND P(n) of the multiphoton squeezed state ( |2^ ) plotted for a nonlinearity F(Xi) = Xf. As in the previous 
case, the eigenvalue f3\ — 3, the squeezing parameter r — 0.8, and the coupling 7 = 0.1. 




A complementary explanation of this remarkable feature in terms of areas of overlap between the Wigner functions 
of number states and of multiphoton squeezed states will be postponed to the section devoted to the phase-space 
analysis of the multiphoton squeezed states. 

Of special physical interest, among the expectations and the correlation functions of_any order in the field operators, 



are the mean number of photons n 



and the mean square deviation An 



a' aa' a) — (a' a 



for 



the number operator n = a^a. The expressions of the first two moments of n are rather complicated in the general 
multiphoton squeezed state d26|), but greatly simplify in the multiphoton squeezed state (p7|). In this state, choosing 
for instance with F(Xi) = XfT 4> = (or 8 = n/2) and /?2 = 0, the mean photon number reads 



[(2 + 3 7 2 ) cosh [2r] - 3 7 2 sinh [2r] - 2] 



+ (3 2 (cosh [2r] - sinh [2r]) (l + 6 7 2 + 4 7 2 /3 2 ) 



(30) 



which, reminding that /3 2 = |a| 2 exp (2r), reduces to the known expression for the squeezed coherent states of the 
harmonic oscillator n = sinh 2 r + |a| 2 when 7 — > 0. The mean square deviation reads 



A 2 n 



-4r 



1 + e Sr + 12 7 2 + 48 7 4 + 2e 4r ( 7 2 - l) 



+ 8f3f (1 + (18 + 4e 4r ) 7 2 + 4/3 2 7 2 (4 + e 4r + 42 7 2 + 16/3 2 7 2 ) + 96 7 4 ) ] 



8/3f (1+ (18 + 4e 4r )y + 4^ 7 
In the limit 7 — > the expression ( |31| ) reduces to 

-4r 

A 2 n= — 



1) +8|a| 



2 e 2r 



(31) 



(32) 



which of course coincides with the mean square deviation in the squeezed coherent states of the harmonic oscillator. 
We see from the above expressions (^) ([n]) that the multiphoton squeezed states can exhibit, just as in the linear 
case, either sub-Poissonian or super-Poissonian statistics depending on the assigned values of the parameters. To gain 
a deeper insight of the field statistics it is useful to study the second-order degree of coherence, i.e. the normalized 
second-order correlation function 



ff (2) (0) 



(at 2 a 2 ) 



= 1 



(33) 



whose values allow to distinguish between the different possible statistical regimes. In particular, if g^ 2 '(0) < 1 
the system exhibits sub-Poissonian statistics, while for g^(0) > 1 the statistics is super-Poissonian. In Fig. 3 the 
correlation function g^ 2 '(0) of the squeezed coherent state of the harmonic oscillator is plotted as a function of the 
parameter of squeezing r at f3\ = 3 (obviously here 7 = 0). 
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FIG. 3. The second-order correlation function <3^ 2 '(0) for the multiphoton squeezed state ( p7j ) plotted as a function of the 
squeezing parameter r at 7 = and j3i = 3. It provides the statistical behavior of the harmonic oscillator squeezed coherent 
state. It is exactly Poissonian (g^ 2 \0) = 1) at r = (coherent state). It then goes sub-Poissonian at small enough values of r, 
turns super-Poissonian at about r = 0.8 until it saturates at the value g^ 2 \0) = 3 for large values of r. 



Turning on the additional nonlinear Bogoliubov interaction with the external medium makes the correlation function 
(0) parametrically dependent on the coupling 7. In Fig. 4 we show the second-order coherence plotted as a function 
of r again at j3± = 3 but now with 7 = 0.05 and quadratic nonlinearity F(X\) = X\. 




0.5 L 

FIG. 4. The second-order correlation function g^ 2 \0) for the multiphoton squeezed state ( p7| ) with quadratic nonlinearity 
F(Xi) = X\ plotted as a function of r at /3i = 3 and 7 = 0.05. The plot shows that, at variance with the corresponding g^ 2 \0) 
for the squeezed coherent state of the harmonic oscillator (7 = 0) shown in Fig. 3, the multiphoton squeezed state exhibits 
super-Poissonian statistics at r = 0. 

Fig. 4 shows that the statistical behavior is significantly affected by the presence of the nonlinearity. Besides 
exhibiting super-Poissonian behavior at r = 0, the multiphoton squeezed state with 7 = 0.05 acquires a sub- 
Poissonian statistics in a more restricted range of values of r (roughly from r = 0.2 to r = 0.6) compared to the case 
7 = 0. If we further increase the value of the coupling 7, the multiphoton squeezed state exhibits super-Poissonian 
statistics for all values of r, including r = 0, as shown in Fig. 5. 
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FIG. 5. The second-order correlation function g^ 2 \0) for the multiphoton squeezed state ( p7| ) with quadratic nonlinearity 
F(Xi) = Xf plotted as a function of r at f3\ — 3 and 7 = 0.1. In this case the multiphoton squeezed state exhibits 
super-Poissonian statistics for all values of the squeezing parameter r. 



It is also interesting to study the field statistics in a multiphoton squeezed state by tracking the behavior of the 
second-order coherence as a function of 7 for fixed values of the squeezing parameter r. In Fig. 6 the correlation 
function g^(0) in the state (|27| ) with quadratic nonlinearity F(Xi) = X\ is studied as a function of 7 for (3\ — 3 and 
r = 0.5. 
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FIG. 6. The second-order correlation function g^iO) for the multiphoton squeezed state ( p7| ) with quadratic nonlinearity 
F(Xi) = Xl plotted as a function of 7 at fi\ = 3 and r = 0.5. The crossover between sub- and super-Poissonian behavior 
takes place at about 7 = 0.05. 



The interplay between squeezing and nonlinear deformation is seen most clearly by studying the correlation function 
g^(0) at larger values of the squeezing parameter r, as shown in Fig. 7. For r > 0.8 the multiphoton squeezed state 
(^7]) exhibits super-Poissonian statistics for all values of 7. 



10 



1 . 8 
1 . 7 



1 . 6 




0.2 0.4 0.6 0.8 1 1.2 1.4 

FIG. 7. The second-order correlation function </ 2 '(0) plotted as a function of 7 at y3i = 3 and r = 0.8. for the state (JSt]) 
with quadratic nonlinearity F(Xi) = X\. The multiphoton squeezed states exhibits in this case super-Poissonian statisitics 
for all values of 7. 

VII. PHASE-SPACE ANALYSIS 

The squeezed states of the harmonic oscillator are typical nonclassical states of light. To study the nonclassical 
features of the multiphoton squeezed states and to compare them with those of the squeezed states of the harmonic 
oscillator it is most convenient to perform a phase-space analysis in terms of the Wigner quasiprobability distribution 
W(Xi, X2). It is well known that the Wigner function is positive-defined for the squeezed states of the harmonic 
oscillator. Although well known, it is plotted in Fig. 8 below with specific values of the parameters for later comparison 
(We have relabelled the axes according to "particle" language: X\ = x, X2 = p). 




FIG. 8. The Wigner function of the harmonic oscillator squeezed state, i.e. the multiphoton squeezed state ( |27| ) with 7 = 0, 
at p\ = 3, r = 0.8. The axes have been relabelled according to the "particle" language of phase space: X\ = x, X2 = p. 

The planar section of the same Wigner function is also useful for later comparison, and it is shown in Fig. 9 below. 
We notice that it is centered at X 2 = p = 0, and it is zero outside a small interval [—5, 5] of variation of p. 
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FIG. 9. The planar section of the Wigner function for the harmonic oscillator squeezed state, i.e. the multiphoton squeezed 
state (H) with 7 = 0, at f3i = 3, r = 0.8. 

Positivity of the Wigner function is preserved only by the multiphoton squeezed states with quadratic nonlincarity 
F(Xi) = X\. Even in this instance however, the Wigner function of the state ( p7| ) is displaced, rotated and deformed 
compared to the corresponding distribution with 7 = 0, as shown in Fig. 10 below. 




FIG. 10. The Wigner function of the multiphoton squeezed state ( P7| ) with quadratic nonlinearity F{X\) = X\ , at /3i = 3, 
r = 0.8, and 7 = 0.1. Notice that it is not centered in p = 0, and it is elongated and rotated compared to the purely squeezed 
case. 

More insight can be gained by looking at its planar section. As shown in Fig. 11 the Wigner function of the 
multiphoton squeezed state exhibits an egg-shaped section (an elongated and deformed ellipsis). Therefore the area 
of overlap with the circular section associated to the generic number state of the harmonic oscillator in phase space 
gives rise to a nontrivial interference which is drastically modified with respect to the purely squeezed case originally 
discussed by Schleich, Walls and Wheeler |l3|, |l4|]. This fact is responsible for the behavior of the oscillations in the 
photon number distribution of a multiphoton squeezed state (see Fig. 1). 
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FIG. 11. Planar section of the Wigner function for the multiphoton squeezed state ( p7| ) with F(Xi) = Xf, f3i = 3, r = 0.8, 
and 7 = 0.1. The elliptic section of the Wigner function for the harmonic oscillator squeezed states is deformed into a egg-shaped 
section. 

We see from Fig. 11 that besides the deformation of the elliptic shape, the Wigner function of the multiphoton squeezed 
state with quadratic nonlincarity, although still positive defined, it is translated and rotated compared to the purely 
squeezed case. In particular, it is not any more centered in p = 0. Moreover it is elongated, and it is non zero only 
in the range of negative values [0,-10] of p. The effects discussed above are significantly enhanced if we consider 
multiphoton squeezed states with higher nonlinearities. Starting with the next higher nonlinearity F{X{) = the 
Wigner function of the corresponding multiphoton squeezed state becomes very strongly deformed and acquires also 
negative values in some regions of the phase space. The Wigner function for the multiphoton squeezed state with 
cubic nonlincarity is reported in Fig. 12 below. 
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FIG. 12. The Wigner function for the multiphoton squeezed state ( p7j ) with cubic nonlinearity F(Xi) = Xf, and the same 
values of the parameters f3i = 3, r = 0.8, and 7 = 0.1 as in the quadratic case F(Xi) = X\. Notice that it becomes negative 
along a strip of values of p ranging from —2 to —24 and of x ranging from 1.4 to 2.4. 
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We notice that, compared to the case with quadratic nonlincarity, the Wigner function is more rotated, tending to 
place itself parallel to the p axis. It is also much more elongated, as it is non zero for negative values of p lying in the 
range [0, —25]. The planar section of the Wigner function for F(Xi) = Xf is plotted in Fig. 13 below. 
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FIG. 13. Planar section of the Wigner function for the multiphoton squeezed state (^) with cubic nonlinearity F(Xi) = Xf , 
at values /3i = 3, r = 0.8, and 7 = 0.1. 

We notice that the elliptic section becomes completely deformed into a narrowing "wing" for increasing negative 
values of p. Finally, in the case F(X\) = Xf the Wigner function becomes even more twisted and shows many waves 
and ripples with negative and positive peaks, as shown in Fig. 14 below. 




FIG. 14. The Wigner function for the multiphoton squeezed state ( |27] ) with quartic nonlinearity F(Xi) — Xf, at values 
/?i = 3, r = 0.8, and 7 = 0.1. 

The planar section has the form reported in Fig. 15 below. It shows clearly the structure of the lateral waves and 
ripples of negative and positive peaeks. Besides, the central body is of still narrower "wing" section with further 
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elongation, as the Wigner function is now non zero in the interval [0, —50] of negative values of p. 
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FIG. 15. Planar section of the Wigner function for the multiphoton squeezed state ( p?] ) with quartic nonlinearity F(Xi ) = Xf , 
at values /3i = 3, r = 0.8, and 7 = 0.1. 

VIII. DISCUSSION AND OUTLOOK 

In the present work we have introduced a new class of nonclassical states of light. These new states have been 
obtained via a quadrature-dependent generalization of the linear Bogoliubov transformation hrst introduced by Yuen 
in his seminal work on the two-photon coherent states of light. We have shown that our quadrature-dependent 
Bogoliubov transformation, generated by a nonlinear operator function of the hrst quadrature is canonical under very 
broad constraints on the numerical coefficients of the transformation. The associated eigenstates of the transforma- 
tion dehne a new class of nonclassical states that preserve some basic features of the coherent states of the harmonic 
oscillator such as classical nondispersive motion of the wave packet center. They also exhibit the amplitude squeez- 
ing typical of the squeezed states of the harmonic oscillator. However, the quadrature-dependent transformations 
generated by quadratic and/or higher powers of the hrst quadrature X\ are associated to non quadratic effective 
Hamiltonians that summarize the interaction with media with higher nonlincarities (multiphoton down conversion 
processes) and we thus name them as multiphoton squeezed states. The presence of the quadrature-dependent term 
in the transformation is reflected in the presence of an additional phase factor in the wave function representation of 
these states. The interplay of the additional phase factor with the standard squeezing contributions is responsible for 
the remarkable nonclassical behavior exhibited by the field statistics, in particular by the photon number distribution 
and the second-order correlation function g^(0). The phase-space analysis in terms of the Wigner quasiprobability 
distribution shows that the multiphoton squeezed states are higly nonclassical, as the Wigner function is strongly 
deformed and it acquires negative values for cubic, quartic, and higher nonlincarities. 

We have discussed the transformation generated by a nonlinear operator valued function of the first quadrature 
Xi . It is of course possible to introduce the transformations generated by a nonlinear operator valued function of the 
second quadrature X2 and also by two nonlinear operator valued functions of each quadrature. However the treatment 
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is in these cases slightly more tedious and involved, and will be deferred to a follow-up of the present work. In fact, 
looking at future developments, it will be interesting to explore the possibility, that we have only mentioned here, 
of considering quadrature-dependent transformations generated by operator valued functions more general than the 
simple powers considered in the present work. For instance, considering the transformation generated by a periodic 
function of the first quadrature might allowc to define squeezed states for massive particles interacting with optical 
lattices. 

We should finally mention the possibility, in contexts broader than quantum optics, of exploiting the quadrature- 
dependent Bogoliubov transformation as a technical tool for the introduction of normal modes. In particular, we 
are currently investigating the possibility of applying the quadrature-dependent Bogoliubov transformation to the 
approximate diagonalization of the Hamiltonian of the weakly interacting Bosc gas. The aim is to diagonalize a 
larger portion of the Hamiltonian in the many-body Hilbert space in comparison with the standard diagonalization 
originally performed via the linear Bogoliubov transformation p5[ . 
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